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Several mechanisms have been proposed to generate primordial magnetic fields and it is often 
assumed that magnetic fields are not affected by a sharp drop in ionization rate due to the cosmic 
recombination. We investigate the validity of the assumption by studying the behavior of magnetic 
fields and fiuid motion around recombination. Fluid equations including the effect of recombination 
QQ , are considered for protons, electrons and neutral hydrogens separately, combining Maxwell equa- 

■ tions. We find that the residual ionization rate required for the conservation of magnetic field at 

' cosmological scales is about 10"^'^, which is much smaller than the standard value ~ 10"*. Further 

^SJ , we will show the acceleration of protons and electrons in the process of recombination which con- 

, ^ . pensates the decrease in the carrier of electric current in order to preserve electric current and then 

^ ^ magnetic fields. 

<\ 

CN ■ I. INTRODUCTION 

<N : 

I— 1 1 There have been a lot of studies on cosmological magnetogenesis before recombination. Among them are 
i-C ■ generation from nonlinear cosmological perturbations during radiation-dominated era [1-11], phase transition 
^p-i" [12, 13] and breaking of conformal invariance during inflationary epoch [14-18] (see [19-21] for comprehensive 
Q reviews) . One of their motivations is to give seed fields to explain the origin of galactic magnetic fields through 
the dynamo amplification. Also, primordial magnetic fields are expected to play a important role in structure 
, formation. In this context, it is usually assumed that magnetic fields survive through cosmic recombination. This 
, assumption seems valid because, even though the ionization rate drops to about 10^^ [22] after recombination, 

magnetohydrodynamic approximation is expected to be still valid due to the residual charged particles and 
comoving magnetic fields are then conserved. 
^ ' Because this is a crucial assumption in the studies on magnetic fields in the early universe, its validity 
. is worth investigating in detail and this is the subject of this paper. We will estimate how small the residual 
1^ ' ionization rate could be in order to preserve magnetic fields with cosmological scales. Furthermore, the evolution 
of the relative velocity between protons and electrons will also be studied in detail. The acceleration of the 
relative velocity should accompany the recombination because the conservation of magnetic fields means the 
conservation of electric current and then, noting the decrease in the current carriers, increase in the relative 
velocity. Thus, to see this, we need to follow the evolution of protons and electrons separately. We pursue 
I these problems by solving the equations for three fluids, protons, electrons and neutral hydrogens, including 
the effect of recombination and combining Maxwell equations. Because we are interested in an epoch around 
the recombination, the density fluctuations and velocities of the fluids are so small that linear perturbation 
would be a very good approximation. Further we adopt tight coupling approximation [23] between the three 
^ fluids which is valid when the collision timescale is much smaller than the dynamical timescale. To make the 
5^ problem tractable, the effects of photons will be neglected as is often done in the previous studies on magnetic 
field evolution after recombination [24, 25]. 

In [24], the evolution of primordial magnetic fields after recombination was considered taking ambipolar 
diffusion into account (see also [25-27]). Because neutral hydrogens do not feel Lorentz force, the presence 
of magnetic fields leads to the difference of motion between neutral hydrogens and charged particles. Then 
neutral hydrogens act as effective electric resistivity which results in the diffusion of magnetic fields. This is 
a nonlinear effect and, therefore, will not be included in the present analysis. However, it was shown that 
ambipolar diffusion is not effective in the tight coupling regime under consideration. As we will see, we have 
another type of effective electric resistivity which appears even at the linear order. 

This paper is organized as follows. In section II, we present basic equations of our analysis: continuity equa- 
tions, Euler equations and Maxwell equations. We will neglect general relativistic effects including cosmological 
expansion because we are interested in electromagnetic phenomena deep inside the horizon and gravitation 
would have minor effects there. Instead, we simply adopt Newtonian gravitational potential which should be 
determined by the distribution of cold dark matter, but our final results do not involve gravitational potential 
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so that we do not need to mention cold dark matter and Poisson equation. The basic equations are solved 
in section III using tight coupling approximation. In order to see the acceleration of protons and electrons 
separately, we need to go on to the next to the lowest order. Finally, we will give discussion and summary in 
section IV. 



II. BASIC EQUATIONS 

The continuity equations of number densities are expressed as, 

dtUp + V • [upVp] = -F, (1) 
dtUe + V • {n,v,) = -F, (2) 
dtUH + V • {nHVH) = F, (3) 

where Ua and Va {a = p, e, H) are the number density and velocity of component a, respectively. Here the F 
terms represent the rate of recombination and we see that [up + uh) is conserved with a mean value, 

np + riH w 3 X 10^ cm'^ f 4ir) • (4) 



, 103 , 

With the same F, the momentum continuity equation of each component can be written as, 

dtirUpUpVp) + V • {nipnpVpVp) = enp{E + VpXB) — e^npUe'qdvpe — rUpUpVpSvpH — mpUpV^ — FnipVp, (5) 

dt{m^n^v^) + V • {rrierieV^Ve) = —ene{E + x B) + e'^npTi^rjSvpp — m,f.n^Ve5v^u — meneV$ — Fm,eVe,{6) 
dtimHriHVH) + V • {mHriHVHVH) = rUpUpUpSvpH + m^n^Ve^v^H - rnHnnV^ + F{mpVp + irieVe), (7) 

where TOq, is mass of a particle a, e is electric charge, E and B arc electromagnetic fields, Sv^fj = Va — Vfj is 
the velocity difference between components a and /?, and $ is gravitational potential. The electric resistivity, 
?7, and the collision frequency between the component a and neutral hydrogen, are [28], 

r,= ^^\nK^lO-'^sec[—^\ — (8) 



T3/2 V 10^ / V 10 



Ue = {vf(TeH)nH ~ 10 ^Xh SeC ^ j ' (9) 

= {vf(rpH)nH ~ IQ-^XH sec-1 {^^^ > (10) 

where T is the plasma temperature, In A is the Coulomb logarithm, v^^^ is the thermal velocity of component 
a, j3 = me/mp and OaH is the scattering cross section between a particle a and a neutral hydrogen. Here we 
defined Xp = np/{np + ur) and xh = nH/{np + nn)- Substituting the continuity equations of the number 
densities into the momentum continuity equations, we obtain, 

o 

6 -* -♦ S Tl 7? 

dtVp + {vp ■ V)vp = — {E + VpX B) —5vpe - VpSvpH - V$, (11) 

2 

dtVe + {Ve ■ y)Ve = —{E + VeXB) + ^^SVpe - V^SVeH - V$, (12) 

rUe me 

TTt Tl Tfl Th F 
dtVH + [vh ■ V)vh = ''^—^VpdVpH H ^—^I'eSVeH - V$ H (rUpVp + meVe - mnVH), (13) 

It is convenient to rewrite equations using the following quantities, 

mpUp + niene „ 

nb= , , Snpe = np-ne, (14) 

nip + me 

nipUpVp + rrieneVe ^ ^ 

Vb = ; , dVpe = Vp-Ve, (15) 

nipUp + merie 
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and conversely, 



np = nb 



Vb 



ne = nb- 



■13 



5n. 



(3 5n. 



pe 



1 + /3 (l + /3)2 nb 



Su 



pet 



+ 



13 Sn. 



pe 



l + f3 + rife 



6v. 



pe' 



In terms of the new variables, the electric current can be expressed as, 

p = e{np - Ue) = eSupe, 



j = e{npVp - UeVe) = e 



nb6u 



pe 



Sn 



■pe^b 



^-^SUpedVpe 



rife 



Svp 



(16) 
(17) 

(18) 
(19) 



We write the recombination rate phenomenologically as F = fub where 1// ^ 1 /H is the effective timescale 
of recombination and is assumed to be constant in time here. Actually, / is a function of several quantities 
such as the ionization fraction and depends on time, but the detailed time profiles of Up, and uh are not 
important for our analysis. One can see that this form of F does not violate the charge conservation. In fact, 
we have, from Eqs. (1) and (2), 



dtUb + V • (rifeUfe) = -fnb, 

dtSUpe + V • {SnpeVb + UbSVpe) = 0. 

We solve the above equations of motion and Maxwell equations, 

dtB = -S/ xE, 



(20) 
(21) 



(22) 
(23) 



to investigate the behavior of three fluids and electromagnetic fields around cosmological recombination. To do 

this, we first expand the physical quantities according to cosmological perturbation theory. At the zeroth order, 
the universe is homogeneous and isotropic so that fluids do not have peculiar velocities and electromagnetic 
quantities vanish. Only the number densities have nonvanishing values and, because wc arc neglecting the 
cosmological expansion, they are constant in time in the absence of the recombination term. At the first order, 
the universe becomes inhomogeneous and anisotropic, and we have peculiar velocities and electromagnetic 
quantities as well as the number densities which arc functions of space and time. Actually, for the allowed 
amplitude of primordial magnetic fields, < 10~^ Gauss (comoving), by CMB observations [29-32], magnetic 
fields can be treated as perturbations. We consider the equations of motion, (11) - (13), and Maxwell equations, 
(22) and (23), up to first order in cosmological perturbations. Wc sec that in this approximation the advection 
terms and the Lorentz force terms of equations of motion are neglected. It should also be noted that the electric 
current is approximated as j w eubSvpe- Then the equations of motion and Maxwell equations reduce to 



dtVp 



dtVe 



and 



e ^ e^UbT] 



rup 



rup 



m 



p 



m. 



■p 



Vp SVpe - Vp5vbH - V$ 



E H 6Vpe - VeSVeH - V$ 

rUe rUe 



dtVH = 



/3 



-7(z/p - Ue)SVpe + 7 



J^p + f3l^e 



1 + /3 



SVbH - V$, 



dtE = V X _B — eubSvpe, 
dtB = -VxE, 



(24) 

(25) 
(26) 



(27) 
(28) 
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respectively. Here 7 = Up/nn and evolves from 00 to 7end ^ 10~^. It should be noted that because we are 
neglecting the second and higher order terms, the coefficients of the velocity terms are zeroth-order quantities. 
Prom Eqs. (20) and (3), the evolution equations for the zeroth-order number densities are, 



rib ' riH 

Further we rewrite the equations of motion in terms of velocity differences: 



dtvb 



dtSVpe 



dtSvbH 



(Vp - VejOVpe , , a OVbH - V$, 



(l+/3)e 



E 



me. 



(l + /3)2 



SVpe - {Vp - l'e)SvbH, 
SvbH- 



-^;^(l+7)+7/ 



(29) 

(30) 
(31) 
(32) 



Equations for the velocity differences, (31) and (32), and Maxwell equations, (27) and (28), consist our basic 
equations to investigate the behavior of magnetic fields and other quantities around cosmological recombination. 



III. TIGHT COUPLING APPROXIMATION 



We solve the equations presented in the previous section using tight coupling approximation [9, 10, 23]. This 
approximation is based on the fact that the different components are tightly coupled through frequent scattering 
so that their velocities are almost the same, that is, 

v^=v + v^^^ +v2'^ +■■■ , (33) 

where v is the common velocity and the others are the deviation from it. The deviation is suppressed by a small 
parameter which is typically the ratio of scattering timescale to dynamical timescale Tdyn- Here the dynamical 
timescalc can be considered as that of acoustic oscillation, which is comparable or several orders less than 
Hubble time, Tdyn « 2 x 10^'^ sec((l + 2;)/10'^)~'^/^. Thus tight coupling approximation is valid because 

the timescale of Coulomb scattering, 

and those of scatterings involving neutral hydrogens, Eqs. (9) and (10), are much smaller than the dynamical 
timescale. Eq. (33) shows that the velocity differences start from the first order, 

5v„, = + 5^1^ +■■■ = (.1^) - «f + (e' - «T0 + • • • • 

Then we see that, in the equations of motion (31) and (32), the time derivative term is higher order than the 
collision term. 

Electromagnetic fields also start from the first order because the exact tight coupling of protons and electrons 
do not allow their existence. 

E = E^'^ + E^"^ +■■■ , (35) 
B = B(^) + B^"^ +■■■ . (36) 

We will solve Eqs. (27), (28), (31) and (32) up to second order in tight coupling approximation. Here it should 
be noted that this expansion is independent of the expansion with respect to cosmological perturbations argued 
in the previous section. For further details of tight coupling approximation, see [9, 10]. 
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A. tight coupling approximation I 



The evolution equations of velocity differences and Maxwell equations at the lowest order in tight coupling 
approximation are, 



(l + /3)e 



me 



l + l-i 



5^1) - (y^ - ^,)5^5<3 = 0, 



:(1+7)K-Z.e)54i)- 



aB(^) = - V X 



^^(1+7) + 7/ 



(37) 

(38) 

(39) 
(40) 



Prom now on, we neglect the 7/ term in Eq. (38) because it is much smaller than v^. We can express and 
^(^) by bv^) using Eqs. (37) and (38), 



/3 V, 



p 



The latter equation is an effective Ohm's law and the effective electric resistivity is. 



rjeS = 11 + 



VpVe 



VpVeTC 



e^ribrj Vp + f3ve Vp + five 



VeTC ~ 2 X 10"S"^- 



(41) 
(42) 

(43) 
(44) 



We can see that the electric resistivity is mainly contributed from Coulomb collision in the early phase, 7^1, 
and collision between charged particles and neutral hydrogens can dominate in the late phase if 7end < 10~^. 
In the realistic case with 7end ~ 10~^, the contribution from neutral hydrogen is at most comparable to that 
from Coulomb collision. 

To solve Eq. (39), we compare the time derivative of electric field, and the velocity difference between 
protons and electrons, bv^} , using Eq. (42), 



' ^-77l-(^d;n + /)^«l 



\enbbiv>P 



dyn 



(45) 



where we replaced the time derivative with Tdyn- Thus, E'^^^ in Eq. (39) is higher order and should appear at 
the next order. Then, we have approximately. 



eub^^^} = V X 



and accordingly, from Eqs. (41) and (42), 



l + (jvp + enb 
E(^) =r?effV X B(^). 



V X B^^K 



(46) 

(47) 
(48) 



Thus wc have expressed the velocity differences and electric field in terms of magnetic field. The evolution 
equation for magnetic field is given by Eq. (40), 



(49) 



This is a diffusion equation and the comoving diffusion length, below which magnetic fields diffuse due to the 
electric resistivity in a Hubble time, is 



Adiff = - \/W « 5 X lO-Vl + 2 X 10-47-^ PC. 



(50) 



6 



The recombination increases the effective electric resistivity by a factor (1 + 2x10 ^7 ^) and then the diffusion 
length by a factor i/l + 2 x 10~'*7~^. Therefore, the evolution of magnetic fields at cosmological scales are 
not practically affected by recombination with such a large residual ionization rate 7end ~ 10~^. According to 
[5, 6, 8-10], magnetic fields are generated by cosmological perturbations with scales from lOOAU to the horizon 
scale. To preserve magnetic fields of these scales, it is enough to have a residual ionization rate as low as 
~ 10-1". 

On the other hand, the diffusion time for a fixed coherence length is, 



H 



5 X 10 



20 



feff 



VlMpc 



j (1 + 2 X 10- V^) ^• 



(51) 



From Eqs. (46), (47) and (48), we see that the velocity differences and electromagnetic fields decrease with 
the timescale (51) which is much larger than cosmological timescale. Because magnetic fields are practically 
constant in time, electric currents which support them must also be constant in time. This can be seen in 
Eq. (46), and this also means that the velocity difference between protons and electrons increases by a factor 
7end ^ IQ-^. It would be instructive to see how protons and electrons are accelerated during the process of 
recombination. However, equations of motion, Eqs. (24), (25) and (26), at the lowest order are. 



dfV 



dtv = 



dfV 



nip 



1 



nie 1 + (3 



(52) 
(53) 
(54) 



all of which lead to. 



dtv = -V$. 



(55) 



Thus, at this order, all components behave as a single fluid and the interactions among the three components 
vanish effectively. This is exactly expected from lowest-order tight coupling approximation. Therefore, we must 
go on to the next order where deviations among the three components begin to appear. 

B. tight coupling approximation II 

The evolution equations of velocity differences and Maxwell equations at the next to the lowest order in tight 
coupling approximation are. 



d,SviO = 



(l + /3)e 



eW - 



(1 + P)e^nbr] _^ pv^ + 



1 + /? 



- i^p - ^^e)5w^H^ 



,(1 +7)K - -e)Sv</P - ^^^(1 +7)<542^ 



= - V X E^"l 



(56) 

(57) 

(58) 
(59) 



Agam we solve Eqs. (56) and (57) for SvIh and E^"^ 



_ P ^p-^e xAII) I /3(t^p - Ve) 



eUbrjeSOV^e 



Svi"^ + 



(l+/3)e 



/3(fp - Vef 



(zVp + /Jz/e)2(l+7) 



dtdvi 



pe 



(60) 
(61) 
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We do not need to solve further for 5vp"' to obtain effective equations of motion for 
substituting Eqs. (60) and (61) into Eqs. (24), (25) and (26), we have 

/3 



and In fact, 



rup 

P 
1 + /? 

TOe 
1 



+ 



1 + 7 

1 V, 



V 



■ l3Ve 



pe 



l+jVp + l3v, 



erib 



(62) 



1 + 
1 

13 7 



1 + 

1 + 13 



me 






Vp - Vf. 


1 + 7 


Vp + (3Ve 




Vp - Z/e 



1 + 7 + 



(63) 



1 + /? 1 + 7 + 



7 



l + /31 + 7i'p + (3ve eub 
In the last equality of each equation, we used, 



(64) 



(65) 



because rjcsk"^ H < f. We sec that purely second order terms have canceled and effective accelerations are 
expressed just by B'^^K The net acceleration of protons is contributed from electric fields and collision with 
neutral hydrogens and they correspond to the first and second term in (• • • ) in Eq. (62), respectively. The 
former is dominant in the early phase of recombination while the latter is dominant in the late phase. The 
acceleration is proportional to the recombination rate / and inversely proportional to the baryon density Ub- 
Thus the acceleration of protons increases substantially as recombination proceeds. The situation is the same 
for electrons except that the first term contributed from electric fields is always dominant. Thus we can see that 
there is no acceleration without recombination, as it should be, and we could confirm that charged particles 
are accelerated to compensate for the decrease in their number densities and maintain the electric currents. It 
should be noted that the direction of acceleration is always the opposite between protons and electrons and the 
compensation is mainly realized by the acceleration of electrons. 

For completeness, we will derive the evolution equation of magnetic fields at this order. Prom Eq. (58), 



eubSv^p'P = V X B(") - 

= V X - r?,2ffV'V X B(^) 



+ /3i/, 



^(l + 7)/VxB(^). 



(66) 



Substituting this into Eq. (59) gives the evolution equation for magnetic fields, 



dtB^"') = -enbTjesV x Sv^JJ^ - 

= VefiV'B^") - 
1 I'pVe 



rUe 



(l+/3)e 
1 

(1 + 7)r?eff 



1 + 



P{Vp - Vef 



Ills 



1 



(z.p+/3l/e)2(l + 7) 

{yp + l3v.f{\ + ^) 

f3{Vp - Ve) 



Picking up the dominant term, we finally have, 

dtB(")=r,,sV'B(") + 



{lyp + Pveni + j) 



dtV X 



/ 



(67) 



(68) 
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IV. DISCUSSION AND SUMMARY 



In our analysis above, we neglected photons and let us mention a couple of potentially important effects. One 
is effective electric resistivity due to Thomson scattering between charged particles and photons. The relative 
importance of the effective resistivity 77x11 and the normal resistivity rj is expressed as [10], 



T] {l + pye-^ribri \ 10-^ ) 



b/2 

(69) 



where ctt is Thomson cross section and p^^ is the photon energy density. Thus, we can see that the effective 
resistivity due to Thomson scattering play a minor role for a reasonable residual ionization rate. On the other 
hand, photons also affect the dynamics of charged particles. In particular. Ohm's law will have a contribution 
from Thomson scattering other than the effective resistivity. The analysis of this term requires a four fluid 
treatment, which is beyond the scope of this paper. 

In this paper, we followed the evolution of protons, electrons, neutral hydrogens and magnetic fields around 
recombination using tight coupling approximation and cosmological perturbation theory. We found that neutral 
hydrogens act as an effective electric resistivity which depends on the ionization rate and that the residual 
ionization rate required to maintain magnetic fields at cosmological scales is about 10"^", which is much smaller 
than the standard value 10~^. Further, we showed that charged particles, especially electrons, are accelerated 
to compensate the decrease in their number density to maintain electric currents. 
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